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ABSTRACT 


The  present  state  of  knowledge  of  the  charged  particle  radiation  environment 
in  synchronous  orbit  is  described.  The  composition,  intensity,  energy  spec¬ 
trum  and  temporal  variations  of  the  radiation  are  listed.  The  data  are  based 
on  measurements  on  board  satellites  in  elliptical  and  in  synchronous  orbits, 
including  the  results  from  the  electron  telescope  on  board  LES-4,  whose 
apogee  is  in  synchronous  orbit.  The  radiation  consists  mainly  of  trapped 
electrons  and  low  energy  protons,  but  solar  flare  protons  also  make  a  sig¬ 
nificant  contribution  The  time-averaged  data  are  given  in  a  form  useful  for 
predicting  the  long-term  radiation  environment  to  which  a  satellite  in  syn¬ 
chronous  orbit  would  be  exposed. 
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CHARGED  PARTICLE  RADIATION  ENVIRONMENT 
IN  SYNCHRONOUS  ORBIT 


I.  INTRODUCTION 


This  report  describes  the  eurrent  state  of  knowledge  of  the  natural  eharged  partiele  radiation 
environment  to  which  a  satellite  traveling  in  an  earth  synchronous  orbit  would  be  exposed.  In 


order  to  include  satellites  in  quasi-synehronous  orbits,  the  region  considered  lies  between 
18,000  and  19,363  nmi.  This  may  be  expressed  in  other  units  as  shown  below. 


Altitude 

nmi  stat.  mi  km 


Radial  Distance 
km 


Earth  Radii 


(,V 


18,000  20,718  33,358 

19,363  22,237  35,884 


39,736 

42,262 


6.23 
6.  63 


A  satellite  in  synchronous  orbit  is  subject  to  the  following  types  of  natural  charged  particle 
radiation  trapped  particles,  mainly  electrons  and  low  energy  protons,  solar  flare  particles, 
mainly  protons;  and  cosmic  radiation.  Of  these,  the  trapped  partieles  are  by  far  the  most  nu¬ 
merous,  but  solar  flare  protons  make  an  important  contribution,  particularly  during  a  sunspot 
maxim  uni . 

For  a  complete  characterization  of  the  radiation  environment,  the  particle  flux  intensity, 
direction  and  energy  spectrum  should  be  measured  as  a  function  of  time,  since  all  these  param¬ 
eters  are  subject  to  both  short  and  long  term  variations.  Ideally,  all  measurements  should  be 
made  from  satellites  in  the  synchronous  orbit. 

The  ATS  1  satellite,  whieh  was  launched  into  the  synchronous  orbit  on  6  December  1966, 
carries  several  experiments  on  board  to  measure  both  the  proton  and  electron  distribution 
(see  Table  I).  Results  giving  a  fairly  comprehensive  description  of  the  radiation  environment 
in  the  synchronous  orbit  are  now  becoming  available.*  * ^  In  the  absence  of  such  data  it  was 
necessary  to  eollate  individual  measurements  from  satellites  in  elliptical  orbits  as  they  tra¬ 
versed  the  synchronous  orbit.  The  comparison  of  these  results  was  complicated  by  temporal 

14 

variations  of  the  fluxes  and  by  differences  in  energy  response  of  the  various  instruments. 

In  the  absenee  of  sufficient  data  from  satellites  in  the  synchronous  orbit,  measurements 
obtained  from  a  satellite  spending  a  large  fraction  of  its  time  in  Synchronous  orbit  is  of  con¬ 
siderable  interest.  FES-4  was  launched  on  21  December  1965  into  an  elliptical  orbit  with  an 
apogee  near  synchronous  orbit  and  a  perigee  of  200  km.  For  30  percent  of  its  time  it  was  lo¬ 
cated  at  R  6.  An  electron  telescope  on  board  the  satellite  monitored  electrons  in  three 
integral  energy  ranges  >150  keV,  >850  keV  and  >2.7  MeV.  It  also  included  a  proton  channel  in 
the  range  1.8  to  50  MeV.  The  data  obtained  from  FES-4  are  compared  with  other  experimental 
results  in  this  report. 


1 


TABLE  1 

RADIATION  DETECTORS  ON  BOARD  ATS-1  SATELLITE  (Ref.  1) 

Orgonizer 

Detector 

Particles 

Energy  Ronge 

Aerospace  Corporation 

Omnidirectional 

Spectrometer 

Electrons 

Protons 

>300  keV,  >450  keV,  >1.9  MeV 
5-21  MeV,  21  -70  MeV 

University  of  Minnesota 

Pratan  Electron 
Spectrometer 

Protons  and 
Electrons 

50  —  150  keV,  150  -500  keV, 
50-1  000  keV 

University  of  Colifornia, 
San  Diego 

Omnidirectionol 

Porticle  Detector 

Protons 

Electrons 

>12  MeV,  >20  MeV 
>0.5  MeV,  >1  MeV 

Bell  Telephone 
Loboratories 

Multi-element  Porticle 
Detecting  Telescope 

Protons 

Alpho  Porticles 
Electrons 

0.7-100  MeV 

1,8-85  MeV 
>1  MeV 

Rice  University 

Suprathermol  Ion 
Detector 

Positive  Ions 
Electrons 

0.25  -50  eV 
>3  keV 

II.  TRAPPED  ELECTRONS 

A.  Vette  and  Lucero's  Data  for  Electrons  at  Synchronous  Altitudes 
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Vette  and  Lucero  have  processed  data  from  Explorers  6,  12  and  14,  Imp  A,  OGO-A,  and 
ERS-17  to  obtain  the  most  comprehensive  information  about  the  electron  distribution  in  the  syn¬ 
chronous  orbit  covering  the  energy  range  between  40  keV  and  2  MeV  and  the  time  range  from 
near  solar  maximum  through  solar  minimum. 

They  found  that  the  time-averaged  omnidirectional  flux  for  energy  greater  than  E  may  be 
expressed  by  an  exponential  function 

JOE,  B/Bq,  *)  -  CA( 4> )  ( B/Rq) "b  exp  l  E/EJ  e/cm2  •  sec  (1) 

for  energies  greater  than  lOkeV,  where 

R  =  the  geomagnetic  field,*  *  ' 

=  the  B-coordinate  at  the  magnetic  equator 

=  local  time  in  hours 

7  7 

C  =  9  x  10  at  solar  minimum  and  4.5  x  10  at  solar  maximum 

with  an  estimated  error  of  2 

b  =  0.625 

E  =  0.215  MeV 
o 

and  A  and  N  are  functions  of  the  local  time.  A  more  useful  presentation  of  the  integral  spectrum 
averaged  over  local  time  at  solar  minimum  is  shown  in  Fig.  1.  At  energies  above  0.5  MeV  this 
spectrum  may  be  expressed  in  the  form 


2 


sec 


(2) 


JOE,  B/Bq)  =  5.04  X  107(B/Bo)'°'625  exp  [-E/0.215]  e/cm2 

At  lower  energies  the  flux  increases  above  the  value  given  by  Eq.  (2).  Table  II  shows  the  average 
omnidirectional  flux  at  the  magnetic  equator  =  1)  at  solar  minimum.  The  figures  are 

believed  to  be  accurate  within  a  factor  of  two.  The  probability  that  the  flux  will  exceed  the  av¬ 
erage  flux  has  also  been  computed,  but  for  missions  with  a  duration  in  excess  of  several  weeks, 

19 

the  average  flux  should  be  used.  W.  D,  Brown  has  compared  Vette  and  Lucero's  model  to 
minimum  and  maximum  flux  rates  measured  by  Aerospace  Corporation  on  ATS-1  in  synchronous 
orbit,  as  shown  in  Fig,  2.  The  combined  Explorer  12  and  14  spectra  are  also  shown.  Consider¬ 
ing  the  factor  of  2  uncertainty  in  Vette  and  Lucero’s  time-averaged  environment,  the  curves  are 
in  good  agreement. 


TABLE  II 

TRAPPED  ELECTRON  FLUX 

Energy  Range 
(MeV) 

Average  Omnidirectional  Flux  (e/cm^  •  sec) 

B/Ba  =  1,  Sunspot  Minimum  (Ref.  15) 

Minimum  Shielding 
ta  Stop  Electrons 
af  Energy  E] 

(in  Al,  Ref.  18) 
(mg/cm2) 

L  =  6.  6,  19, 363  nmi 

L  =  6.0,  1 

7,300  nmi 

Integral 
(E  >E,) 

Differential 

(E  |  <  E  <  E^) 

Integral 

(E  >  E  j) 

Differential 

(E,<E<E2) 

E. 

E2 

0.01 

0.02 

6.51  X  107 

7.1  X  106 

5. 69  X  107 

2.  2  X  io6 

0.25 

0.02 

0.05 

5.  80  X  107 

1.  15X  107 

5. 46  X  107 

6.  1  X  106 

0.85 

0.05 

0.  10 

4.65X  107 

1.  15X  107 

4.85X  107 

8.9  X  106 

4.3 

0.  10 

0.20 

3.  50  X  107 

1 . 40  X  107 

3.96X  107 

1. 30  X  107 

14.  1 

0.20 

0.40 

2.  10  X  107 

1.31  X  107 

2.66X  107 

1.47X  107 

41. 8 

0.40 

0.60 

7.  94  X  106 

4.87X  106 

1.  19X  107 

6.5  X  106 

158 

0.60 

0.80 

3.07X  106 

1.88X  106 

5.37X  106 

2.96X  106 

210 

0.80 

1.00 

1.  19  X  106 

7.2  X  105 

2.41  X  106 

1.33X  106 

309 

1.00 

1.20 

4.66X  105 

2.  84  X  105 

1.08X  106 

5.9  X  105 

424 

1.20 

1.40 

1 . 82  X  105 

1.11  x  io5 

4.87X  105 

2.68X  105 

530 

1.40 

1.60 

7.  14  X  104 

4.34X  104 

2.  19  X  105 

1.21  X  105 

636 

1.60 

1.80 

2.  80  X  104 

1. 70  X  104 

9.83X  104 

5.41  X  104 

742 

1.80 

2.00 

1.  10X  io4 

6. 7  X  103 

4. 42  X  104 

2.44X  104 

848 

2.00 

2.50 

4.  32  X  104 

3.  90  X  103 

1 . 98  X  104 

1. 71  X  104 

954 

2.56 

3.00 

4.  18  X  102 

3.  77  X  102 

2. 69  X  103 

2. 33  X  103 

1220 

3.00 

3.50 

4.06X  10 

3-66  X  10 

3.63X  102 

3.  14  X  102 

1490 

3.50 

4.00 

3.95 

3.57 

4.92X  10 

4.25X  10 

1750 

4.00 

5.00 

3.  84  X  10_1 

3.  80  X  IO’’ 

6.66 

6.54 

2010 

-3 

-1 

5.00 

3.64X  10 

1. 22  X  10 

2540 
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B.  LES-4  Data 


Energetic  particle  fluxes  were  measured  by  an  experiment  aboard  LES-4.  This  instrument 
consists  of  several  surface  barrier  transmission  detectors  stacked  and  operated  with  coincidence 
techniques  in  order  to  identify  particle  nature  and  energy.  The  response  of  the  various  channels 
to  monoenergetic  electrons  was  measured  as  a  function  of  energy.  The  criterion  for  selection 
of  threshold  energy  was  minimum  change  in  effective  efficiency  with  spectral  characteristic 
energy.  Only  3  channels  will  be  considered  at  this  point  with  selected  thresholds  at  150  keV, 

850  keV,  and  2.75  MeV. 

The  omnidirectional  flux  has  been  obtained  by  suitable  averaging  of  directional  measurements. 
The  low  energy  channel  has  been  corrected  for  bremsstrahlung  while  coincidence  techniques 
eliminate  the  problem  in  the  upper  channels.  The  >2.75  MeV  channel  has  been  cor  rected  for 
galactic  cosmic  rays  penetrating  the  shielding. 

Points  consisting  of  5-minute  averages  were  found  for  every  observed  passage  of  LES-4 
through  the  L  =  6.6  region  from  launch  through  the  next  747  days.  The  average  over  these  56 
points  is  plotted  in  Fig.  3.  Both  local  time  and  magnetic  storm  effects  should  be  smoothed 
because  of  random  sampling.  B/Bq  ranged  from  1  3  to  5  in  these  samples. 

Also  plotted  are  spectra  from  the  points  in  which  the  850  keV  channel  had  both  its  maximum 
and  minimum  count  at  L  =  6.6. 

Comparison  of  these  results  to  those  shown  in  Fig.  2  indicate  quite  good  agreement. 

C.  Variation  of  Electron  Flux  with  L 

It  is  desirable  to  ascertain  the  variation  of  the  electron  flux  and  of  the  energy  spectrum  with 

geocentric  equatorial  distance  Rq  for  satellites  in  a  quasi-synchronous  orbit. 

If  the  earth’s  geomagnetic  field  were  that  of  a  perfectly  centered  dipole,  the  magnetic  shell 

L  in  the  L,  B  coordinate  system  developed  by  Mcllwain^’^  should  be  equal  to  R^  measured  in 

units  of  earth  radii.  In  the  vicinity  of  the  earth’s  synchronous  orbit,  the  geomagnetic  field  is 

disturbed  by  magnetic  storms  and  the  steady  state  solar  wind  and  undergoes  both  diurnal  and 
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other  irregular  temporal  variations.  ’  The  location  of  the  L-shell  is  then  no  longer  fixed 
but  varies  with  changes  in  the  earth’s  magnetic  field.  The  experimental  data  should  therefore 
be  more  properly  plotted  against  Rq>  but  in  this  report  the  parameter  L  will  be  retained  follow¬ 
ing  the  sources  from  which  the  data  have  been  obtained. 

Figure  4  shows  the  average  and  the  maximum  omnidirectional  flux  at  L  =  6.0,  L  =  6.6  and 
L  =  7.0.  The  results  are  based  on  calculations  by  J.  I.  Vette  on  OGO-A  spectrometer  data  in 

the  energy  range  from  50  to  690  keV  obtained  by  Winckler  and  Pfitzer  of  the  University  of 
22  —  26 

Minnesota.  The  average  values  at  L  =  6.6  and  L  =  7.0  show  an  enhancement  of  the  low  en¬ 

ergy  electron  flux  in  excess  of  the  exponential  spectrum.  The  energy  spectrum  softens  with  in- 
27 

creasing  L.  The  average  omnidirectional  flux  at  L  =  6.0  may  be  approximated  to  an  exponential 
spectrum  of  the  form 

JOE)  =  5.91  x  107  e"E//0-25  e/cm2  •  sec  .  (3) 

The  constant  5.91  X  107  has  been  normalized  for  the  average  flux  at  the  magnetic  equator  at 
solar  minimum  by  comparison  with  the  constant  of  Eq.  (2).  Values  of  J  for  L  =  6  are  also  shown 
in  Table  I. 


4 


Figure  5  is  a  comparison  of  the  variation  of  electron  flux  with  L  from  data  obtained  by 
Explorer  14  (Ref.  15)  and  LES-4.  The  Exp.  14  data  have  been  corrected  to  local  noon  and  repre¬ 
sent  the  median  values  over  the  lifetimes  of  the  data.  The  LES-4  data  are  simple  means  over 
the  data  at  different  L  values.  The  B/B  and  local  time  variation  is  similar  for  different  L' s. 

Since  e.  very  definite  variation  with  L  occurs  it  should  be  pointed  out  that  a  geostationary 
satellite  would  be  expected  to  experience  an  average  flux  dependent  on  the  longitude  at  which  the 
satellite  is  fixed.  The  variation  in  L  with  longitude  is  shown  in  Fig.  6. 


D.  Temporal  Variations 

The  flux  level  of  the  outer  zone  electrons  varies  in  time  by  about  two  orders  of  magnitude. 

As  an  example.  Fig.  7  indicates  the  values  of  the  fluxes  of  >150keV  and  >850  keV  electrons  as 

sampled  by  LES-4  at  L  =  6  over  a  period  of  200  days.  The  detectors  on  board  the  ATS-t  satel- 

1—13  28 

lite  have  provided  detailed  information  on  these  variations. 


1.  Diurnal  or  Local  Time  Variations 


The  distortion  of  the  geomagnetic  field  in  the  synchronous  orbit  by  the  solar  wind  produces 

29 

an  average  diurnal  variation  in  the  electron  intensity.  The  particles  at  midnight  are  on  particle 
shells  which  lie  outside  those  sampled  at  noon  and  hence  have  a  lower  particle  intensity.  The 
magnitude  of  the  variation  is  a  function  of  the  electron  energy  with  the  largest  decreases  in  flux 
on  the  night  side  appearing  for  the  high  energy  electrons.  The  electron  spectrum  is  softer  at 
night  than  during  the  day 

During  quiet  and  moderately  disturbed  conditions  the  maxima  and  minima  in  the  range  from 

4  5 

300  to  1900  keV  are  displayed  about  an  hour  ahead  of  local  noon  and  midnight  ’  On  the  other 

28 

hand,  Freeman  and  Maguire  found  that  thermal  plasma  ions  show  maximum  fluxes  in  the 
midnight  to  dawn  sector. 

During  quiet  times  there  is  a  smooth  diurnal  variation  which  becomes  less  smooth  as  mag¬ 
netic  activity  increases.  During  more  disturbed  times,  the  satellite  leaves  the  region  of  stable 
trapping  on  the  night  side  and  enters  a  part  of  the  magnetosphere  which  is  inhabited  by  rapidly 
fluctuating  electron  populations.^  On  some  occasions,  typically  following  a  magnetic  storm,  no 
diurnal  variation  is  found;  instead  the  flux  increases  in  a  stepwise  fashion  throughout  the  post 
storm  recovery  period. 

Typical  diurnal  changes  during  quiet  or  moderately  disturbed  periods  are  as  follows; 


E  ^  300  keV 
E  ^  450  keV 
E  1050  ke\ 
E  £  1900  keV 


Maximu m  Variation 
3 

3.5 

9 

20 


Minimum  Variation 
1.5 
1.5 
3 
3 


The  variations  are  not  correlated  to  K 
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Before  the  results  from  the  ATS-1  satellite  became  available  Yette  and  Lucero  adopted 

an  iterative  procedure  to  obtain  from  all  their  data  the  composite  local  time  variation  shown  in 

29  30 

Fig.  8.  Measurements  of  the  local  time  variation  reported  earlier  by  Frank  ’  are  included. 
The  figure  shows  that  the  local  time  variation  increases  with  energy  and  peaks  around  9  o'clock 
with  a  minimum  at  midnight.  The  range  of  variations  of  the  data  obtained  from  ATS-1  appears 
to  be  somewhat  smaller  than  that  of  the  data  from  elliptical  orbits.  This  is  mainly  due  to 
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differences  in  data  selection  and  analysis.5  The  data  obtained  during  very  disturbed  days  have 
been  omitted  in  the  analysis  of  ATS-1  results. 

2.  Short  Term  Changes 

Some  small  scale  periodic  fluctuations  of  electron  fluxes,  with  a  period  of  100  —  1000  seconds 

depending  on  the  energy,  are  associated  with  the  longitudinal  drift  of  bunches  of  electrons  that 

1 3 

have  been  injected  into  the  magnetosphere. 

During  magnetic  storms  the  electrons  >45keV  may  also  undergo  quasi-regular  or  irregular 
modulations  by  one  order  of  magnitude  with  a  typical  period  of  six  minutes.  * 

3.  Variations  Due  to  Magnetic  Storms 

At  L  =  4  there  is  a  strong  correlation  of  electron  flux  with  magnetic  storms.  Each  large 

flux  increase  is  accompanied  by  an  increase  of  spikiness  of  K  .  An  initial  decrease  in  flux  at 

the  onset  of  the  storm  is  followed  by  a  rapid  rise  during  the  main  phase  of  the  storm  by  about 

two  orders  of  magnitude  with  a  rise  time  of  10  minutes  or  longer  depending  on  the  energy.  A 
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slow  exponential  decay  then  sets  in  which  continues  until  the  onset  of  the  next  storm. 

In  the  synchronous  orbit  the  electron  flux  is  also  affected  by  conditions  on  the  sun,  but  the 
variations  are  far  more  irregular  and  do  not  correlate  well  with  magnetic  storms.  This  is 

9 

believed  to  be  due  to  physical  motion  of  the  magnetospheric  boundary.  Electrons  in  the  energy 
range  50  —  150  keV  show  large  increases  of  intensity  that  persist  for  30  —  90  minutes  and  repeat 
with  a  period  of  2  —  4  hours.  The  largest  intensity  occurs  in  the  midnight  to  dawn  sector. 

Energetic  electrons  show  a  greater  response  than  the  low  energy  ones.  Electrons  >1.9  MeV 
decrease  by  one  order  of  magnitude  disappearing  at  local  midnight.  This  is  followed  by  a  gradual 
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recovery  which  takes  place  more  slowly  than  for  lowr  energy  electrons. 

Occasionally,  the  synchronous  altitude  moves  right  outside  the  magnetosphere,  usually 
within  three  hours  of  local  midnight.  Subsequently  discrete  clouds  of  energetic  particles  are 
detected  which  execute  repeated  longitudinal  drifts  about  the  earth,  as  well  as  bursts  of  very 
low  energy,  E  <  50  keV,  highly  directional  positive  ions.7 

4.  7-Day  and  27-Day  Variations 

3 

A  near  periodic  7-day  recurrence  of  flux  depressions  has  been  observed  similar  to  the  low 

3 

altitude  variations  related  by  Williams  to  changes  in  the  interplanetary  field  direction. 

Magnetic  activity,  and  hence  the  electron  flux  in  the  synchronous  orbit,  has  a  tendency  to 

recur  after  the  lapse  of  one  or  more  intervals  of  27  days,  the  approximate  period  of  rotation 

3  3 

of  the  sun  relative  to  the  earth  moving  along  its  orbit.  The  27 -day  variations  can  persist  for 
long  time  periods  like  a  year  but  are  not  always  present. 

5.  Variation  with  Solar  Cycle 

15 

Vette  and  Lucero  found  no  evidence  for  an  increase  in  the  trapped  radiation  flux  levels  at 

L  =  6.6  as  solar  maximum  is  approached.  Their  data  showed  a  decrease  by  a  factor  2,  which  is 

within  the  intercomparison  accuracy  of  the  various  satellite  measurements.  This  is  an  agreement 

34 

with  O’Brien’s  observation  that  the  range  of  intensities  of  the  outer  zone  electrons  in  1963  was 
the  same  as  in  1959,  even  though  the  solar  activity  strongly  influences  the  outer  zone  intensities 
via  magnetic  storms.  Presumably,  the  principal  effect  of  the  solar  maximum  vs  solar  minimum 
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will  simply  be  an  increase  in  frequency  (and  to  a  lesser  extent  in  magnitude)  of  changes  of  the 
outer  zone. 

35 

Frank  and  Van  Allen  found  that  the  inner  edge  of  the  outer  radiation  belt  moved  outward 
during  the  declining  phase  of  solar  activity  cycle  19.  The  changes  produced  by  this  effect  in  the 
spectral  distribution  of  the  electron  flux  at  L  =  6.6  appear  to  be  within  the  margin  of  error  of 
different  experimental  observations. 

The  LES  4  instrument  sampling  in  the  synchronous  region  shows  no  indication  of  a  dependence 
of  flux  levels  on  solar  activity.  Data  examined  include  some  800  days  of  operation,  wherein  the 
smoothed  sunspot  number  changed  from  20  (near  its  minimum  value)  to  100  (above  the  maximum 
of  the  mean  of  sunspot  cycles  1  through  19). 

E.  Pitch  Angle  Distribution 

The  pitch  angle  distribution  determines  the  direction  from  which  the  trapped  charged  par¬ 
ticles  will  strike  a  satellite  in  orbit.  The  pitch  angle  is  a  function  of  the  particle  energy  and 
undergoes  temporal  variations  like  all  the  other  parameters. 

Pitch  angle  distributions  for  electrons  in  the  synchronous  orbit  have  been  reported  by 
Lezniak  et  al.  ^  The  distribution  at  local  noon  during  a  disturbed  period  is  shown  in  Fig.  9.  The 

low  energy  distribution  is  highly  anisotropic.  Similar  measurements  are  believed  to  have  been 

23  —  26 

carried  out  by  Pfitzer  and  Winckler  on  OGO  A.  The  measurement  may  be  made  with  a 

single  directional  detector  if  the  satellite  spin  axis  is  perpendicular  to  the  magnetic  field.  Other¬ 
wise  a  series  of  directional  detectors  are  required  covering  the  quadrant  from  0  to  90°.  The 
direction  of  the  spin  axis  must  be  known  with  respect  to  the  magnetic  field.  In  the  synchronous 

orbit  the  pitch  angle  distribution  cannot  be  derived  from  the  omnidirectional  intensity  along  a 

3  6 

line  of  force,  since  the  adiabatic  invariance  conditions  are  not  obeyed. 

F.  Low  Energy  Electrons 

37-39 

Frank  measured  the  intensity  of  electrons  in  the  range  from  200  cY  to  50  kcY  near  the 

magnetic  equator  with  an  array  of  electrostatic  analyzers  on  OGO-3.  The  electron  energy  den- 

—8  3  37 

sities  remained  fairly  constant  at  3  x  10  erg/cm  over  L  ^  4.0  —  7.5  corresponding  to  an 

9  2 

integral  electron  flux  >200  e\  on  the  order  of  10  e  cm  •  sec.  This  is  in  good  agreement  with 

9  2 

an  integral  elec  tron  flux  %6l0eV  of  3  X  10  e/cm  •  sec  at  a  geocentric  distance  of  8  Rq  (Fig.  6, 
Ref.  24).  At  higher  energies  the  flux  varies  considerably  with  radial  distance  (Fig.  1,  Ref.  38). 

The  differential  energy  spectrum  decreases  monotonically  with  energy  at  1.  3.9  (Fig.  10,  Ref.  39), 

whereas  at  10. 1  R  it  shows  a  broad  peak  at  2  keV  (Fig.  8,  Ref.  38).  No  spectral  data  have  been 
published  closer  to  the  synchronous  orbit. 

40  4  1 

Earlier  measurements  of  low  energy  electrons  at  L  -  6.6  were  made  by  Freeman  *  on 

Explorer  12  using  cadmium  sulfide  total  energy  flux  detectors  and  also  by  Russian  experimenters 
42  4  3 

on  Zond  2  and  Mars  1  employing  charged  particle  traps.  Their  results  were  more  qualitative 

in  nature,  but  are  in  broad  agreement  with  Frank’s  data. 

44  -46 

Serbu  and  Maier  have  measured  the  electron  distribution  from  0  to  100  e\  near  the 

geomagnetic  equatorial  plane  by  means  of  retarding  potential  analyzers  on  IMP-1  and  -2.  The 
electrons  and  positive  ions  below  4eV  constitute  a  neutral  thermal  plasma  with  a  Maxwell- 
Boltzmann  energy  distribution.  In  the  synchronous  orbit  the  electron  density  is  in  the  25  to 
80  e/cm  range  and  the  electron  temperature  lies  between  1  and  2eV.  The  electron  flux  in 
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9  2 

the  0  <  E  <  5  eV  energy  range  was  found  to  be  3  X  10  e/cm  •  sec  on  IMP-1.  There  is  also  a 

higher  energy  component  whose  flux  in  the  5  <  E  <  10  eV  energy  range  is  2  x  10  e/cm  •  sec. 

47 

Slysh  has  derived  the  plasma  density  from  antenna  impedance  measurements  on  Zond-2. 
He  obtained  a  density  of  100  e/cm  at  6R  and  a  radial  dependence 

N  =  1.3  X  105  R  '4 
e  o 

valid  from  4  to  7  R^  in  good  agreement  with  the  results  of  Serbu  and  Maier. 

in.  TRAPPED  PROTONS 

A.  Intensity  and  Energy  Distribution 

The  positively  charged  trapped  particles  in  the  outer  radiation  zone  consist  predominantly 

48  49 

of  low  energy  protons.  Van  Allen  and  Krimigis  *  have  recently  found  trapped  alpha  particles 
with  a  maximum  intensity  at  L  ~  3.1.  The  r^tio  of  the  differential  spectral  intensity  of  alpha 
particles  to  that  of  protons  is  of  the  order  of  0.05  at  2  MeV,  and  the  ratio  of  the  integral  direc¬ 
tional  intensities  (E  >  0.52  MeV/nucleon)  is  2.3  x  10 

Proton  measurements  in  the  synchronous  orbit  were  made  by  Davis  and  Williamson^0  ^ 

53  54 

on  Explorers  12  and  14,  by  Vernov,  et  al.,  on  Cosmos  41  and  by  Krimigis  and  Armstrong 

on  Mariner  4. 

Davis  and  Williamson  measured  protons  in  the  energy  range  from  100  keV  to  4.5  MeV  with 

a  directional  scintillation  counter.  Vernov,  et  ad.,  used  semiconductor  counters  sensitive  to 

protons  in  two  ranges,  0.4  MeV  <  E  <  7  MeV  and  3  MeV  <  E  <  8  MeV.  Krimigis  and  Armstrong 

also  measured  protons  in  two  energy  ranges,  0.50  4  E^  ^11  MeV  and  0.88  E^  <  4  MeV  on  the 

geomagnetic  equator.  The  experimental  data  of  Davis  and  Williamson  have  been  processed  by 
55 

J.  H.  King  of  Aerospace  Corporation,  together  with  data  from  proton  measurements  in  lower 
orbits  which  cover  the  energy  range  from  400  keV  to  4  MeV.  From  these  data  a  model  of  the 
proton  environment  has  been  constructed  which  gives  the  omnidirectional  flux  for  energy  greater 
than  E  MeV, 

JOE,  B,  L)  =  F(B.  L)  exp  [{0.4  -  E)/E  (B,  U]  p/cm2  •  sec  .  (4) 

Figure  10  shows  the  energy  distribution  of  the  protons  at  L  =  6.6  as  a  function  of  B/Bq.  The 
values  of  the  distribution  function  F  and  the  spectral  function  Eq  at  the  magnetic  equator  are: 


L  =6.6 

L  =  6.0 

F(Bo,  L) 

4.00  x  105p/cm^  •  sec 

1.6  x  lO^p/cm^  •  sec 

Eo(Bo'  L) 

0.111  MeV 

0.112  MeV 

Calculations  by  J.  I.  Vette  on  experimental  data  by  Davis  and  Williamson  extend  the  proton  dis¬ 
tribution  down  to  100  keV  as  shown  in  Fig.  11.  There  is  a  significant  enhancement  of  the  flux  at 
energies  below  250  keV  over  the  exponential  distribution.  Table  III  shows  the  average  omni¬ 
directional  flux  at  the  magnetic  equator  obtained  from  Eq.  (4)  for  energies  above  400  keV  and 
from  Fig.  11  for  energies  between  100  and  400  keV.  The  accuracy  of  the  proton  model  environ¬ 
ment  is  estimated  to  be  within  a  factor  of  two  for  L  ^  4.5.  Larger  amplitude  temporal  variations 
and  a  relative  shortage  of  data  decrease  the  accuracy  at  synchronous  orbit. 
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TABLE  III 

TRAPPED  PROTON  FLUX 

Energy  Range 
(MeV) 

2 

Omnidirectional  Flux  (p/cm  •  sec) 

B  =  B 

a 

Minimum  Shielding 
ta  Stop  Protons  af 
Energy  E^ 

(in  Cu,  Ref.  56) 
(mg/cm2) 

L  =  6. 6,  1 9,363  nmi 

L  =  6.0, 

1  7, 300  nmi 

1 ntegral 

(E  >  E,) 

Differential 
(E,  <E<E2) 

Integral 

(E>E1) 

Differential 
(E,  <  E  <  E2) 

E1 

E2 

0.  10 

0.25 

1 . 2  X  107 

1.0X  io7 

0.  16 

0.25 

0.50 

1. 7X  106 

1.5X  106 

0.51 

0.50 

0.  75 

1 . 6  X  105 

1. 5X  105 

6.  4  X  105 

5.  9  X  105 

1.45 

0.  75 

1.00 

1.7X  104 

1.4  X  104 

6.  8  X  104 

6.  1  X  104 

2.8 

1.00 

1.25 

1. 7X  103 

1. 5X  103 

6.8  X  103 

6.  1  X  103 

4.5 

1.25 

1.50 

1 . 8  X  102 

1 . 5  X  102 

7.  1  X  102 

6.  4  X  102 

6.  75 

1.50 

1. 75 

1.8X  10 

1.6X  10 

7.1  X  10 

6.  3  X  10 

8.  1 

1. 75 

2.00 

1.9 

1.7 

7.7 

6.9 

11.5 

2.00 

2.25 

1 . 9  x  10_1 

1. 7X  10"' 

7.  7  X  10"' 

6.  9  X  IO'1 

14.5 

-2 

-2 

-2 

-2 

2.25 

2.50 

2. OX  10 

1.8X  10 

8. OX  10 

7.  1  X  10 

17.5 

-3 

-3 

-3 

-3 

2-50 

2.75 

2.  2  X  10 

2. 0  X  10 

8.  7X  10 

7.8X  10 

21 

-4 

-4 

-4 

-4 

2.75 

3.00 

2.2X  10 

2. OX  10 

8.  7X  10 

7.9X  10 

24.5 

-5 

-5 

3.00 

2.  1  X  10 

8.  4  X  10 

28 

B.  Temporal  Variations 


During  magnetic  storms,  changes  occur  in  the  proton  intensity  with  apparent  temporary 

enhancement  of  proton  flux  during  the  main  phase.  At  L-values  between  5  and  8.  large  time 

variations  do  occur  for  protons  with  energies  greater  than  1  MeV,  but  protons  near  too  keV  are 
2  7 

relatively  more  stable. 

Davis  and  Williamson52  (Fig.  12)  show  the  effect  of  a  solar  event  on  5  October  1962,  High 
energy  protons  970  ke\  are  sufficiently  enhanced  to  produce  a  separate  peak  at  I  5.8.  The 
520-keV  protons  are  also  enhanced,  whereas  the  >141-keV  protons  are  only  slightly  disturbed. 
Six  days  later  the  enhancements  iiave  decayed.  The  proton  buildup  is  coincident  with  high  latitude 
negative  bays. 

Data  obtained  trom  Kxplorer  14  (Fig.  13)  show  that  above  1.7  MeV  the  proton  intensities 

fluctuate  by  at  least  one  order  of  magnitude  at  L  6.  whereas  low  energy  proton  intensities  are 

52 

much  more  stable.  Davis  and  Williamson  compared  the  proton  intensities  measured  in  19t>5 
with  those  measured  in  1962.  They  found  that  t lie  intensity  of  protons  above  1.7  !Yle\  throughout 
the  region  from  L  =  2.5  to  5  was  larger  by  a  factor  of  two  to  three  in  1965.  The  enhancement 
extended  to  energies  down  to  500  ke\  .  whereas  100  to  300  ke\  protons  were  greatly  depleted. 

Also,  the  1962  Davis-Williamson  data  (">0.5  Mev)  are  in  apparent  disagreement  with  the  data 
obtained  in  1964  on  Mariner  4.  The  latter  is  about  two  orders  of  magnitude  smaller  than  the 
former  for  L  ^6.0  55 
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28 

Freeman  and  Maguire  have  made  a  detailed  study  of  the  temporal  variations  of  low  energy 
protons  in  the  synchronous  orbit.  During  periods  of  moderate  magnetic  activity,  a  high  particle 
flux  is  seen  near  local  midnight  and  the  particle  intensities  after  midnight  are  higher  than  those 
before  midnight.  During  periods  of  high  magnetic  activity,  the  particle  fluxes  undergo  large 
local  time  variations.  The  enhanced  particle  flux  distribution  extends  over  a  larger  fraction  of 
the  night-side  portion  of  the  synchronous  orbit  and  the  pre  and  postmidnight  asymmetry  is  some¬ 
what  reduced. 

C.  Pitch  Angle  Distribution 

Davis'^  has  measured  the  equatorial  pitch  angle  {EPA')  distribution  at  L  =  6  for  proton 
energies  of  97  and  495  keV  (Fig.  14).  The  distribution  is  fairly  isotropic  with  a  maximum  at  90° 
and  is  subject  to  temporal  variations.  Data  obtained  on  different  passes  show  even  larger  dis¬ 
crepancies,  particularly  for  protons  with  energies  greater  than  495  keV. 

D.  Low  Energy  Protons 

37  39 

Frank  ’  measured  the  low  energy  proton  distribution  near  the  magnetic  equatorial  plane 

in  the  range  200  eV  <  E  <  50  keV  with  an  array  of  sensitive  electrostatic  analyzers  mounted  on 

OGO-3.  The  omnidirectional  proton  intensity  at  L  =  6  in  the  energy  range  1.8  <  E  <  48  keV  was 
8  2 

found  to  be  2{±1)  x  10  p/cm  ■  sec.  The  proton  energy  density  at  L  =  6.5  in  the  range 

—  8  3 

190  eV  <  E  <  48  keV  is  3.9  X  10  erg/cm  at  a  geomagnetic  latitude  of  22°.  This  is  about  4 

5 1 

times  the  energy  density  observed  by  Davis  for  protons  of  energy  ^97  keV. 

The  above  data  represent  quiet-time  conditions.  During  the  main  phase  of  a  geomagnetic 

storm  the  proton  intensities  at  L  6  to  6.5  in  the  energy  range  200  eV  ^  E  50  keV  increase  by 

a  factor  2  to  2.5.  The  proton  angular  distributions  during  the  storm  were  found  to  be  almost 

isotropic  with  maximum  directional  intensities  at  an  equatorial  pitch  angle  of  90°  (Ref.  23). 

The  proton  differential  energy  spectrum  at  L  =  6.0  and  6.6  is  shown  in  Fig.  15.  The  data 

39 

were  obtained  from  Figs.  5,  6,  7  and  1 1  of  Frank's  paper  and  include  higher  energy  values 

calculated  by  Frank  from  the  integral  energy  spectrum  at  L  -  6.0.  EPA  =  65°  published  by 
5  1 

Davis.  The  following  features  are  worth  noting 

A  very  good  fit  of  Frank  and  Davis'  data  at  higher  energies  despite 
an  interval  of  almost  3  years  between  the  measurements. 

At  about  1  keV  the  proton  flux  decreases  below  background  response 
of  the  instrumentation. 

Below  1  keV  the  proton  flux  decreases  as  L  increases,  whereas 
above  10  keV  the  opposite  is  true. 

46 

Serbu  and  Maier  have  measured  the  distribution  of  positive  ions  from  0  to  45  eV  by  means  of 
a  retarding  potential  analyzer  on  IMP-2.  The  positive  ions  have  a  Maxwell-Boltzmann  energy 
distribution  with  a  temperature  of  about  3  to  5  eV  at  a  geocentric  distance  from  5  to  7  earth  radii. 
The  density  of  positive  ions  in  this  region  declined  rather  sharply  from  300/cm'*  at  Rg  =  5  to 
20/ cm'*  at  R ^  =  7. 

The  Rice  University  ion  detector  aboard  the  ATS-1  satellite  has  measured  a  highly  directional 
flux  of  positive  ions  moving  with  a  bulk  flow  velocity  of  30  km/sec  toward  the  sun.  During  one 
event  the  alpha/proton  density  ratio  was  approximately  3%.  the  proton  density  6  protons/cm3 
and  the  temperature  1  eV. 
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Ion  mass  spectrometer  measurements  on  OGO-A'^7  have  yielded  a  ratio  of  100  to  1  in  the 
He4  to  H+  concentrations  over  a  range  of  altitudes  from  1500  to  30,000  km. 

IV.  SOLAR  FLARE  PARTICLES 
A.  Solar  Flares 

Solar  flares  occur  in  a  limited  region  of  the  chromosphere,  usually  in  the  vicinity  of  an 
active  sunspot  group.  The  phenomena  produced  by  a  solar  flare  as  observed  from  the  vicinity 
of  the  earth  may  be  divided  into  simultaneous  and  delayed  effects  as  shown  in  Table  iv  ^8,59 


TABLE  IV 

SOLAR  FLARE  EFFECTS 

Effect 

Energy 

Type 

Travel  Time 

Simul  toneous 

8.  3  minutes 

Visible  Light 

Rodio  Woves 

Ultraviolet  ond  X-roys 

Deloyed 

Salar  Cosmic  Ray  Porticles 

10  MeV  -  >1  GeV 

Mainly  protons 

D  irect 

<1  hour 

Indirect 

20  minutes  ta 

20  hours 

Magnetic  Storm  Particles 

20  ta  50  hours 

Energetic  Storm  Porticles 

100  keV-  10  MeV 

Moinly  protons 

Solar  Flare  Plosmo 

10  eV  -  100  keV 

Ians  and  electrons 

The  importance  of  a  solar  flare  is  defined  in  terms  of  the  area  occupied  by  the  flare  at  the 
time  of  maximum  brightness The  most  important  solar  flares  (importance  classes  3  and  3  +  ) 
occur  very  rarely,  last  for  periods  up  to  8  hours  and  contribute  most  of  the  solar  particles. 
Conversely,  flares  of  importance  1  to  1—  (subflares)  occur  very  frequently,  have  durations  of 
the  order  of  minutes  and  usually  do  not  contribute  any  significant  particle  fluxes. 

Most  of  the  largest  solar  flare  events  produce  direct  cosmic  ray  particles,  i.e..  particles 
from  a  fairly  well-defined  source  in  the  general  vicinity  of  the  sun.  The  time  of  arrival  of  the 
particles  is  commensurate  with  rectilinear  travel  from  the  sun  and  the  intensity  increases  to  its 
maximum  value  in  3  to  20  minutes.  The  indirect  radiation  is  isotropic  in  nature  and  possesses 

a  steeper  energy  spectrum;  also,  the  travel  time,  the  rise  to  peak  intensity  and  the  decay  take 

59 

far  longer.  In  many  events  only  indirect  radiation  is  observed.  The  energetic  storm  particles 
do  not  arrive  in  appreciable  numbers  until  about  a  day  after  a  solar  flare  and  reach  a  maximum 
at  the  time  of  arrival  of  the  magnetic  storm  at  the  earth  which  coincides  with  the  arrival  of  the 
solar  flare  plasma. 
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B,  Composition 


1.  Positively  Charged  Cosmic  Hay  Particles 

The  composition  of  solar  cosmic  rays  has  been  reviewed  by  Biswas  and  Fichtel The 
primary  radiation  outside  the  atmosphere  consists  predominantly  of  protons,  but  also  contains 
alpha  particles  and  heavy  nuclei.  The  energy  spectrum  of  the  flare  particles  varies  with  each 
solar  event  and  also  over  the  course  of  a  given  event.  The  low  energy  spectrum  is  further  in¬ 
fluenced  by  the  geomagnetic  cutoff  at  the  measuring  point  and  by  interaction  with  the  atmosphere 

if  the  measurement  is  made  on  the  ground  or  at  low  altitudes. 

(y  2 

The  Solar  Proton  Manual  states  that,  at  least  for  proton  energies  >30  MeV,  the  average 
differential  energy  spectrum  over  the  course  of  an  event  may  be  expressed  in  the  form 


d.J 

dE 


KE 


3.5  <  n  <5 
o 


(5) 


b3 

Freier  and  Webber  have  shown  that  a  better  representation,  not  only  for  protons  but  also  for 
alpha  particles  and  medium  nuclei  with  A/Z  -  2,  is  given  by  an  exponential  rigidity  spectrum 


=  Jq(  A.  Z,  t)  exp  (-  P/Po(A,  Z.  t)  ] 


(6) 


where  P  is  the  magnetic  rigidity,  or  momentum  per  unit  charge  =  m^vc/Zc,  which  is  related 
to  the  kinetic  energy  E  (in  electron  volts)  by 
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2  2  2  2 
P  Z  =  K  +  2m  c  E 
o 


(7) 


Zc  =  particle  charge 
A  =  atomic  number 
t  =  time 


m  -  rest  mass 
o 


v  -  velocity 
c  -  velocity  of  light 


For  protons,  m^e  =  938  MeV 


P  lies  in  the  range  from  50  to  200  MV;  varies  from  102  to  10^  p/cm2  •  sec  for  solar 

protons.^  This  spectrum  applies  down  to  rigidities  of  0.2GV.^  Both  the  constant  K  and  the 

characteristic  rigidity  P  vary  in  a  systematic  manner  during  the  course  of  an  event.  In  most 

events  there  is  an  exponential  rise  to  a  maximum  intensity  1  r  followed  by  a  slower  exponential 

m  ax  ^  y 

decay.  The  integrated  intensity  above  a  given  energy  may  then  be  calculated  from 

J(>R)  =  £  >max(>K)  exP<  t/tR)  dt  +  [  Imax(>E)  ^P^V  dt  (8) 

where  t  0  is  the  time  of  the  peak  intensity,  t^  and  t^  are  characteristic  rise  and  decay  times 
which  are  energy  dependent  and  vary  from  several  hours  to  several  days  in  different  events. 

The  quantity  t  depends  on  the  location  of  the  flare  on  the  sun;  t.  is  influenced  more  by  general 

I* 

interplanetary  conditions.  The  ratio  of  t^/t^  varies  from  1.6  to  6  with  an  average  of  4. 
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The  integrated  intensities  also  exhibit  exponential  rigidity  spectra  of  the  form 

FOE)  =  G  exp  [-P(E)/P  ]  (9) 

The  characteristic  integrated  rigidity  PQ  is  close  to,  but  slightly  smaller  than,  that  of  the  spec¬ 
trum  at  peak  intensity.  The  solar  protons  produced  by  the  flare  of  23  February  1956  exhibited 
a  P  of  260  MV.  With  this  exception  P  varies  from  50  to  200  MV  Webber^  plotted  the  number 
of  events  having  a  P^  greater  than  a  given  rigidity  P  as  a  function  of  P  for  all  major  events 

from  1956  to  1962.  This  is  shown  in  Fig.  16.  The  linear  relationship  indicates  that  P  varies 

69  ° 

randomly  with  a  constant  probability  density  and  a  mean  value  of  97  MV,  There  is  no  correla¬ 
tion  between  Pq  and  the  total  integrated  flux  or  the  flare  location  on  the  sun’s  surface. 

Alpha  particles  and  heavier  nuclei  possess  approximately  the  same  rigidity  spectrum  as 
the  protons.  The  ratio  of  protons  to  alphas  above  a  fixed  rigidity  lies  most  commonly  between 
one  and  two  although  events  with  ratio  greater  than  50  have  been  observed.  \\  obber^  has  ob¬ 
served  a  correlation  between  the  proton/alpha  particle  ratio  and  the  characteristic  rigidity  P 

as  shown  in  Fig.  17.  The  ratio  increases  with  P  .  A  similar  correlation  between  this  ratio 
**  o 

70 

and  the  cosmic  particle  energies  was  noted  by  Sakurai.  All  the  heavy  nuclei  appear  to  have 
the  same  spectral  shape  and  relative  abundances  shown  in  Table  V  taken  from  the  review  by 
Biswas  and  Fichtel^*  which  also  shows  the  corresponding  figures  for  galactic  cosmic  rays. 
Among  the  heavier  nuclei  carbon  6,  nitrogen  7,  oxygen  8,  and  neon  10  are  the  most  abundant. 


TABLE  V 
(Ref. 61) 

RELATIVE  ABUNDANCES  OF  NUCLEI  NORMALIZED 

TO  BASE  OF  1.0  FOR  OXYGEN 

ENERGY  RANGE:  42  —  135  MeV/NUCLEON 

E  lement 

Solar  Cosmic  Rays 

Galactic  Cosmic  Rays 
(10  ta  30%  uncertainty) 

?He 

107  ±  14 

48 

3Li 

- 

0.3 

4Be-5B 

<0.02 

0.8 

6c 

0.59  ±0.07 

1.8 

7n 

0.  19  ±0.04 

<0.8 

8o 

1.0 

1.0 

9f 

<0.03 

•SO.  1 

10Ne 

0.  13  ±0.02 

0.30 

11  No 

- 

0.  19 

'2Mg 

0.043  ±0.011 

0.32 

,3ai 

- 

0.06 

14Si 

0.033  ±0.011 

0.  12 

15P-21Sc 

0. 057  ±0.017 

0.  13 

22Ti-28Ni 

^0.02 

0.28 

_  3 

Very  small  quantities  of  deuterons  and  tritons  have  also  been  observed,  probably  10  or  less 
of  the  proton  abundance. 

71  2 

M.  V.  K.  Apparao,  et  al.,  have  detected  150  solar  neutrons/m  •  sec  •  sr  with  energy  be¬ 
tween  20  and  160  MeV.  The  main  source  of  neutrons  is  the  neutron  albedo  generated  by  knock-on 
collisions  of  solar  flare  and  cosmic  ray  energetic  particles  with  atmospheric  particles  over  the 

magnetic  poles.  About  10  percent  of  all  neutrons  generated  in  this  way  escape  into  space.  This 

3  2 

resulted  in  albedo  neutron  fluxes  in  synchronous  orbit  during  1967  of  8.6  X  10  n/cm  •  day  in 

4/2  19 

the  absence  of  solar  activity  and  of  1  72  x  10  n/cm  ■  day  during  solar  flares.  These  values 
are  comparable  to  the  solar  neutron  intensity  measured  by  Apparao,  et  al. 

2.  Magnetic  Storm  Particles 

The  magnetic  storm  particles  may  be  divided  into  two  major  rigidity  groups  with  an  exponen¬ 
tial  rigidity  spectrum  [Eq.  (6)]  and  with  the  following  parameters:^ 

Characteristic  Rigidity  P  Flux  J 

_ _ _ J  o  _ o 

5  7  2 

Energetic  Storm  Particles  5  to  30  MV  10  to  10  /c.nri  •  sec 

9  2 

Solar  Flare  Plasma  0.5  to  2  MV  10  /cm  sec 

The  energetic  storm  particles  are  thus  far  more  numerous  than  the  solar  cosmic  rays,  but 
possess  far  lower  energies.  At  10  Me\  the  two  particle  populations  have  approximately  equal 
intensities.  The  energetic  storm  particles  consist  mainly  of  protons,  but  their  exact  composi¬ 
tion  is  yet  to  be  determined.  Considerable  numbers  of  electrons  are  also  associated  with  mag¬ 
netic  storms  (see  below). 

Energetic  storm  particles  with  the  same  rigidity  spectrum  but  with  very  much  lower  intensity 
3  2 

(J  from  10  to  10  particles/cm  •  sec  )  are  also  observed  during  small  magnetic  storms  which 

°  72 

recur  over  many  27 -day  cycles  of  the  solar  rotation  period.  A  major  solar  flare  starts  each 

7  3  74 

series  of  recurrent  events.  Bartels  ’  postulated  that  these  storms  are  produced  by  plasma 
streams  which  are  continuously  emitted  from  long-lived  regions  of  the  sun  called  M -regions/ 5 
Long-lived  solar  plasma  streams  with  a  27 -day  structure  were  observed  by  Mariner  II.  The 
recurrent  storm  particles  are  produced  over  a  period  of  1  to  2  days  with  the  intensity  increasing 
and  then  decreasing  more  or  less  symmetrically  in  time  about  the  time  of  peak  intensity.  They 
consist  mainly  of  protons,  but  helium  nuclei  are  also  present  /  7 

3.  Electrons 

Energetic  electrons  are  very  rare,  much  less  than  one  percent  of  the  protons  at  1  Me\  . 

Relatively  large  fluxes  of  low  energy  electrons  (<1  MeV)  have  been  observed.  Peak  fluxes  of 

up  to  40,000  e/cm7  •  sec  of  energy  greater  than  40keV  have  been  measured  in  interplanetary 

7  8 

space  following  solar  chromospheric  flares.  (  The  electrons  belong  to  two  classes:  prompt 
electrons  which  arrive  within  an  hour  of  the  flare  and  delayed  electrons  which  appear  at  earth 
about  20  to  40  hours  following  a  large  solar  flare  and  are  present  from  8  to  24  hours.  The  prop¬ 
agation  of  energetic  solar  flare  electrons  is  usually  highly  anisotropic  and  takes  place  in  cones 
of  propagation;  outside  the  boundaries  of  these  cones,  the  electron  flux  is  much  reduced.  The 
cones  of  the  prompt  electrons  possess  opening  angles  of  either  30’  (simple  events)  or  90  (com¬ 
plex  events).  The  buildup  time  of  the  prompt  electron  is  less  than  one  hour;  the  delay  time  is 
of  the  order  of  a  day.  The  delayed  electrons  0.5  MeV  possess  comparable  peak  intensities. 
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4.  X-Rays 
79 

DeJager  has  suggested  one  classification  of  solar  x-rays  by  energy: 

Soft  x-rays  \  >  1  A  or  E  <  10  keV 

Hard  x-rays  \  <  1  A  or  E  >  10  keV 

and  another  by  the  mechanism  of  origin: 

Quasi-thermal  x-rays  emitted  by  a  hot,  stationary  plasma  (up  to  10/qK) 

Nonthermal  x-ray  bursts  due  to  the  acceleration  of  energetic  jets  of 
charged  particles 

Gamma  rays  of  nuclear  origin. 

No  solar  gamma  rays  have  been  detected.  An  upper  limit  for  the  flux  of  gamma  rays  above 

50  MeV  is  10  ^  photons/cm  •  sec/*0  The  nonthermal  bursts  consist  of  short-lived  bursts  of 
4  6 

hard  x-rays  (10  to  10  eV)  with  durations  ranging  from  a  few  seconds  to  a  few  minutes  and  with 

an  energy  flux  up  to  2  X  10  ^ergs/cm^  •  sec.^  ^ °  The  spectral  shape  appears  to  be  independent 

of  intensity  and  decreases  rapidly  from  20  to  150  keV. 

Solar  flares  are  also  accompanied  by  changes  in  intensity  and  spectral  composition  of  the 

quasi-thermal  x-rays  emitted  by  the  solar  corona  which  are  summarized  in  Table  VI  taken  from 

8  7 

the  review  article  by  Mandelstam.  Thermal  x-rays  may  be  divided  into  a  quasi-stable  com¬ 
ponent  emitted  by  undisturbed  coronal  regions  and  a  slowly  varying  component  emitted  by  active 

o 

regions.  Below  about  20  A  the  spectrum  is  continuous  and  mainly  determined  by  the  slowly  vary- 

o 

mg  component.  Between  20  and  100  A  the  contributions  of  the  two  components  arc  roughly  equal 
and  a  line  spectrum  produced  by  elements  present  in  the  corona  is  superimposed  on  the  continuous 

CJ 

spectrum.  The  disturbed  corona  subjects  the  intensity  below  20  A  to  considerable  variations  over 

periods  of  a  few  days  and  up  to  several  months.  Disturbed  and  undisturbed  conditions  are  also 

correlated  to  periods  of  maximum  and  minimum  solar  activity.  During  class  2-3  solar  flares, 

the  boundary  of  the  spectrum  is  displaced  down  to  about  1  to  0.  1  A  and  the  flux  below  10  A  rises 

by  an  order  of  magnitude.  The  total  flux  below  100  A  increases  by  a  factor  2.  Figure  18,  taken 

8  8 

from  an  article  by  H.  H.  Malitson,  shows  the  x-ray  spectrum  down  to  1  A  under  different  solar 
conditions. 


TABLE  VI 

SOLAR  X-RAYS  (Ref.  87) 

X 

Quiet  Sun 

Flores 

(A) 

(erg/cm^  •  sec) 

(erg/ern^  .  sec) 

-5  -3 

-2 

,<10 

10  to  10 

10 

-4  -2 

-2 

^20 

10  to  10 

4  X  10 

^100 

1 0-1  to  1 

2 

-6  ^  -5 

^0.1 

10  to  2  X  10 

15 
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Arnoldy.  et  al. ,  established  a  correlation  between  the  integrated  x-ray  flux  with  both  the 
3-  and  10-em  integrated  radio  flux.  The  x-ray  bursts  were  expressed  in  terms  of  absolute 
energy  flux  assuming  an  exponential  x-ray  spectrum  with  an  e-folding  value  of  7  keV.  In  the 
period  from  5  September  1964  to  30  June  1966,  there  were  at  least  3  x-ray  events  with  an  inte- 
grated  flux  between  0.1  and  1  erg/ cm  in  the  energy  range  from  10  to  50  keV. 

C.  Temporal  Variations 

There  is  an  enormous  variation  in  the  frequency  with  which  solar  events  oceur  and  in  the 

integrated  particle  intensity  produced.  In  fact,  the  eight  largest  events  of  solar  cycle  19  have 

produced  about  90  percent  of  all  solar  cosmic  rays  above  10  MeV.  The  annual  integrated  in- 

65  67 

tensities  are  shown  in  Table  VII.  ‘  ''  From  the  limited  data  available,  the  following  conclusions 
may  be  drawn: 

(1)  Almost  the  entire  integrated  intensity  for  the  solar  cycle  is  produced 
during  the  six  years  of  maximum  activity.  The  intensity  accumulated 
during  sunspot  minimum  is  insignificant. 

(2)  The  annual  integrated  intensities  are  lower  during  sunspot  maximum 
(1957  —  1958)  than  during  the  adjacent  periods  of  ascending  and  partic¬ 
ularly  descending  portions  of  the  solar  cycle. 

(3)  The  characteristic  rigidity  PQ  is  smaller  during  sunspot  maximum  than 
during  the  ascending  and  descending  portions  of  the  solar  cycle. 

Satellite  measurements  of  the  integrated  intensities  of  magnetic'  storm  particles  >1  Me\  are 

available  only  since  1961.  Their  mean  characteristic  rigidity  P  is  20  MV.  The  integrated 

°  o  65 

intensities  for  1959  and  I960  have  been  derived  from  polar  glow  measurements  at  3914A. 

Attempts  have  been  made  to  correlate  the  annual  integrated  intensity  of  the  cosmic  ray  par¬ 
ticles  with  the  number  of  class  3+  events,  the  number  of  major  events,  and  with  the  monthly 

smoothed  sunspot  number  averaged  over  the  year.  The  smoothed  sunspot  number  H  for  a  given 
92  ° 

month  is  defined  as 

5 

H0  =  (R-6  +  r+6  +  2  S  'V/24  (,0) 

-5 

where  H .  is  the  observed  monthly  Zurich  relative  sunspot  number  in  months  after  that  month,  in 
which  the  observed  relative  number  was  R  .  The  relative  sunspot  number  is  defined  as 

H  -  k(  lOg  +  f)  (11) 

where  g  is  the  number  of  sunspot  groups,  i.e.,  isolated  clusters  of  sunspots,  f  is  the  total 
number  of  distinct  spots  and  k  is  a  scale  fat  tor  depending  on  the  observer.  For  the  relative 
sunspot  numbers  observed  at  Zurich,  k  =  0.6. 

In  Fig.  19  the  integral  solar  particle  intensity  at  four  different  energy  levels  is  plotted  as  a 

function  of  the  average  annual  smoothed  sunspot  number  S  over  sunspot  cycle  19.*  ^  In  years  in 

which  S  is  less  than  50,  the  total  number  of  solar  particles  above  10  MeV  is  less  than  the  far 

8  2 

more  penetrating  galactic  particle  background  of  ~  10  /cm  •  year. 

The  integrated  proton  intensity  1  MeV  shows  a  fair  degree  of  cor  relation  with  S. 
Approximately. 
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and 


log1Q  1(>1  MeV)  =  9  +  0.02  S 


log10  I(>10MeV)  ■  7  +  0.02  S  .  (12) 

At  higher  energies  there  is  no  correlation.  This  is  not  surprising  since  a  single  flare  may  be 
responsible  for  the  entire  annual  intensity  of  the  high  energy  solar  protons.  Therefore,  the 
correlation  appears  to  be  between  S  and  the  energetic  storm  partic  les  rather  than  the  solar 
protons. 

There  is  also  some  correlation  between  the  number  of  events  N  and  S.  On  the  average 

N  -  0.08  S  (13) 

An  event  is  here  defined  as  one  producing  a  detectable  intensity  increase  of  cosmic  rays  with 
energy  >10  MeV  above  the  galactic  cosmic  ray  background  of  4  partieles/cm^  •  sec  near  earth.00 

The  maximum  sunspot  number  attained  in  a  solar  cycle  has  varied  from  45  in  cycle  6  to 
201.3  in  cycle  19.  The  peak  sunspot  number  in  the  last  four  cycles  has  undergone  a  systematic 
increase.  Predictions  made  before  the  start  of  cycle  20  and  based  on  the  80 -year  term  varia¬ 
tions  in  the  sunspot  maxima  indicated  that  cycle  20  will  exhibit  much  lower  sunspot  numbers 
93-95 

than  cycle  19.  More  accurate  estimates  of  the  magnitude  and  timing  of  the  sunspot  maxi¬ 

mum  of  cycle  20  became  possible  toward  the  end  of  19  07  using  Waldmeier's  relationship  between 
maximum  sunspot  number  R.^  and  ascent  time  t  for  an  even  cycle,  i.e.,  the  time  from  minimum 
to  maximum  sunspot  numbers^ 

log  R  V1  2.69  -  0.17  t  .  (14) 

Webber00  has  predicted  a  maximum  sunspot  number  of  80  to  100  occurring  in  1969.  Figure  20 

97 

shows  the  smoothed  sunspot  numbers  up  to  January  1967.  The  sunspot  minima  have  all  been 
placed  at  the  same  point  along  the  time  axis.  The  smoothed  sunspot  numbers  so  far  observed 
in  cycle  20  indicate  that  the  peak  sunspot  number  will  be  below  that  of  cycles  18  and  19  and 

slightly  above  the  mean  for  cycles  8  through  19.  The  predicted  values  shown  for  cycle  20  have 

98 

been  obtained  by  a  statistical  method  developed  by  McNish  and  Lincoln. 

It  remains  to  be  shown  that  the  correlation  between  sunspot  number  and  annual  proton  in¬ 
tensity  >1  MeV  and  10  MeV  derived  during  cycle  19  remains  valid  during  the  sunspot  maximum 
99 

of  cycle  20.  Ellison  has  compared  cycles  18  and  19.  He  found  that  the  frequency  of  the  large 
flares  determining  the  integrated  solar  proton  intensity  did  not  increase  with  the  sunspot  number 
which  was  made  up  mainly  bv  small  and  medium-sized  spot  groups.  This  suggests  a  rather  loose 
correlation  between  the  sunspot  number  and  the  integrated  proton  intensity.  Also,  in  earlier- 
cycles  showing  lower  sunspot  numbers  the  sunspots  were  concentrated  in  larger  groups  with 
large  penumbral  areas.100  These  are  likely  to  produce  more  intense  flare  particles. 


D.  Geomagnetic  Cutoff 

Tire  earth's  geomagnetic  field  acts  as  a  shield  against  incoming  solar -flare  particles  of  low 

,,  ,  .  ,  r,  ..  ,  ,  r.  ,  ,  10  i  v,  ,  102-  Kb 

energies.  By  applying  the  Stormer  theory  lor  a  magnetic  dipolar  field,  Ray  and  Sauer 

have  derived  a  cutoff  rigidity 


p  =  m.9/l2c;v 

C 


( 15) 


18 


! 


for  a  point  in  space  defined  by  Mcllwain's  magnetic  shell  parameter  L  measured  in  units  of  earth 
16  17 

radii.  ’  Satellite  studies  showed  that  the  effective  geomagnetic  cutoff  rigidity  of  solar  cosmic 
rays  are  significantly  less  than  the  predicted  values  even  during  magnetically  quiet  periods  and 
are  further  reduced  markedly  during  the  main  phase  of  geomagnetic  storms Although 
this  effect  is  variable  from  event  to  event,  Webber^7  has  applied  a  correction  in  the  form  of  an 
average  effective  value  L  ^  in  Eq.  (15)  from  which  the  average  cutoff  rigidity  may  be  obtained. 
Table  VIII  shows  the  average  cutoff  rigidity  and  cutoff  energy  E(  of  solar  protons  and  alpha 
particles  calculated  in  this  manner  for  equatorial  orbits  at  L  =  6  and  L  6.6. 


TABLE  VIII 

AVERAGE  GEOMAGNETIC  CUTOFF 

ENERGIES  (Ref. 67) 

L  =  6.0 

L  =  6.6 

l  r  r 

9.6 

10.5 

eff 

P 

162  MV 

135  MV 

c 

E 

13.9  MeV 

9.  7  MeV 

c 

There  is  considerable  doubt  that  these  cutoff  values  which  have  been  derived  from  observa¬ 
tions  on  the  ground  have  any  validity  in  synchronous  orbit.  In  fact  it  is  known  that  during  the 
peak  of  a  magnetic  storm  the  synchronous  orbit  may  lie  outside  the  earth’s  magnetic  field 
altogether. 

Recent  results  from  the  ATS  -  i  satellite^’^’  *  1  have  shown  that  protons  of  energy  >25  IVIeV 
penetrate  the  synchronous  orbit  at  all  times.  A  variable  shielding  effect  does  exist  in  the  energy 
range  from  6  to  20  Me\  .  There  is  a  diurnal  variation  with  a  maximum  between  1800  and  0400 
local  time,  showing  that  proton  access  is  not  uniform  (see  also  Ref.  108).  It  was  further  observed 
that  protons  of  solar  origin  persist  for  long  periods  of  time  in  the  synchronous  orbit,  so  that  dur¬ 
ing  solar  maximum  interplanetary  space  is  seldom  free  of  them.  Recurrent  streams  of  protons 
are  observed  near  the  middle  of  the  27 -day  solar  rotation  cycle. 

Figure  21  is  a  plot  of  5 -minute  count  rate  av  er  ages  from  LES-4  on  2  September  1966  which 
shows  flare  particles  penetrating  deep  into  the  magnetosphere,  where  one  normally  thinks  only 
of  stable  trapping.  The  electron  channels  >150  keY  and  >850  keV  show  perfectly  normal  behavior. 

At  0648  ITT  the  channel  sensitive  to  protons,  1.8  McY  <  FI  <  50  MeY.  starts  to  rise.  An  optical 
flare  of  importance  3  was  observed  to  start  at  0542.  The  decrease  starting  at  about  0823  is 
presumably  due  to  the  effect  of  the  earth’s  magnetic  field.  Thus,  in  this  case  both  ’’trapped” 
particles  and  particles  which  have  apparently  come  from  the  sun  occupy  the  same  region  of  space. 

It  must  be  concluded  that  there  exists  at  present  no  satisfactory  method  of  calculating  average 
annual  intensities  of  solar  protons  below  10  MeY  in  synchronous  orbit.  Only  an  upper  limit  can  be 
given  assuming  no  geomagnetic  shielding. 

E.  Prediction  of  Solar  Flare  Intensities  in  Synchronous  Orbit 

All  predictions  on  future  solar  flare  particle  intensities  are  based  on  data  obtained  during 
cycle  19,  since  no  meaningful  data  were  available  for  earlier  cycles.  Figure  22  shows  the  integral 
annual  solar  proton  intensity  as  a  function  of  rigidity  for  three  cases  of  interest- 
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(1)  Average  intensity  for  the  years  195  6  —  19  61,  covering  the  years  of 
maximum  solar  activity  of  cycle  19. 

(2)  Intensity  for  1959,  the  highest  intensity  ever  recorded. 

(3)  Intensity  for  1956,  which  includes  the  solar  flare  of  23  February  with 
a  characteristic  integrated  rigidity  of  260  MV,  the  highest  rigidity 
ever  recorded. 

Many  of  the  predictions  of  solar  cosmic  ray  intensities  are  based  on  a  "typical "  solar  flare 
109 

presented  by  Bailey.  This  hypothetical  solar  proton  event  has  been  constructed  from  data 
obtained  in  the  course  of  six  major  events  during  the  last  sunspot  cycle  to  arrive  at  a  time  history 

for  a  typical  large  event.  By  integrating  the  flux  at  different  energies  over  the  time  of  a  solar 

110 

event  Keller  has  obtained  the  integral  and  differential  energy  spectra  associated  with  the 
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"typical  flare"  shown  in  Fig.  23.  A  similar  model  developed  by  Lockheed  is  shown  in  Fig.  24. 
The  two  integral  spectra  are  also  shown  in  Fig.  22  for  comparison  with  the  observed  annual 
spectra.  It  will  be  noted  that  Bailey's  integrated  curve  matches  the  average  annual  spectrum 
from  1956  through  1961  very  closely  both  in  intensity  and  in  the  spectral  distribution.  The  Lock¬ 
heed  spectrum  is  harder  and  closer  to  the  1956  curve. 

By  using  the  regression  coefficients  between  annual  sunspot  number  and  annual  integrated 

solar  particle  intensity  for  cycle  19  [Eq.  (12)],  and  assuming  a  similar  regression  curve  for 
65 

cycle  20,  Webber  arrives  at  the  following  values  for  the  worst  year  of  cycle  20: 

J(>1  MeV)  ^  10  1 1  p/cm2  •  year 

J(>10  MeV)  ^  109  p/cm2  •  year  (16) 

These  values  are  an  order  of  magnitude  less  than  the  maximum  values  for  cycle  19.  As  explained 
in  Sec.  C,  they  should  apply  primarily  to  energetic  storm  particles. 

The  predicted  integrated  solar  flare  proton  spectrum  over  the  period  1967  —  1973  is  shown 
in  Fig.  25.  The  average  curve  is  based  on  the  assumptions  that 

(1)  J(>1  MeV)  =  lO^p/cm2  •  year. 

(2)  Characteristic  rigidity  of  the  energetic  storm  particles  is  the  same 
as  that  of  the  1956  —  1961  average  shown  in  Fig.  15. 

(3)  The  solar  flare  proton  flux  above  30  MeV  is  the  same  as  that  of  the 
1956  —  1961  average. 
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The  average  curve  is  practically  identical  with  a  predicted  curve  by  RCA  also  based  on 
Webber's  data. 

A  maximum  predicted  curve  has  been  constructed  from  the  limiting  conditions: 

(1)  Characteristics  of  the  energetic  storm  particles  are  the  same  as  for 
the  average  curve. 

(2)  The  solar  flare  proton  flux  between  30  and  100  MeV  is  the  same  as 
that  of  the  1959  maximum. 

(3)  At  higher  energies  the  flux  is  the  same  as  that  of  the  195  6  curve, 
with  a  characteristic  rigidity  of  2  60  MV. 

Webber's  two  predicted  points  are  also  shown  on  the  curve.  Geomagnetic  shielding  has  not  been 
considered.  The  alpha  particle  spectra  shown  in  Fig.  26  have  been  derived  from  the  correspond¬ 
ing  proton  spectra  by  assuming  a  proton/alpha  particle  ratio  of  2:1  in  accordance  with  the  value 
for  an  average  value  of  Pq  100  shown  in  Fig.  10.  In  the  maximum  curve,  the  1959  spectrum 
has  been  continued  to  energies  above  100  MeV,  since  in  the  1956  spectrum  the  high  characteristic 
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rigidity  is  coupled  to  a  proton/alpha  particle  ratio  of  50;  1.  The  RCA  spectrum,  which  is  based 
on  a  proton/alpha  particle  ratio  of  10:1,  is  also  shown. 

The  predicted  integral  and  differential  proton  and  alpha  particle  fluxes  are  listed  in  Tables  IX 
and  X  as  a  function  of  energy.  The  average  values  should  be  used  for  long-term  applications  over 
periods  of  several  years.  The  integral  fluxes  at  9.7  and  13.9  MeV,  corresponding  to  maximum 
geomagnetic  shielding  at  L  -  6.0  and  L  =  6.6,  are  also  shown.  The  minimum  shielding  thick¬ 
nesses  have  been  derived  from  the  proton  range-energy  curves  in  copper/*8  assuming  that  alpha 
particles  of  the  same  velocity  have  the  same  range. 


TABLE  IX 

PREDICTED  SOLAR  FLARE  PROTON  FLUX  OVER  THE  PERIOD  1967-1973 


Energy  Range 
(MeV) 

2 

Pratan  Flux  (p/cm  •  year) 

Minimum  Shielding 
to  Stop  Protons 
af  Energy  E] 

(in  Cu,  Ref.  56) 
(mg/ cm2) 

Average 

Maximum 

Integral 
(E  >E,) 

Differential 
(E,  <  E<  E2) 

Integral 
(E  >E,) 

Differential 
(E,  <  E  <  E2) 

E, 

E2 

0.  1 

0.5 

4.2  X  1011 

2.  3  X  lo" 

4. 2  X  101 ' 

2.  3  X  101 ' 

0.  16 

0.5 

1.0 

1.9  X  lo" 

9.  0  X  1010 

1.9  X  1011 

9. OX  lo'° 

1.45 

1.0 

2.0 

1.0  x  lo" 

4. 4  X  lo'° 

1.0  x  io1 1 

4. OX  lo'° 

4.5 

2.0 

3.0 

5.6  X  1010 

2.  3  X  1010 

6.0  X  1010 

1. 7x  1010 

14.5 

„  ,„io 

10 

10 

in 

3.0 

4.0 

3.3  X  10 

1.  1  X  10 

4.3  X  10 

1.0  X  10 

28 

4.0 

5.0 

2.2  X  1010 

8.0  X  109 

3.3  X  1010 

8. OX  109 

46 

5.0 

7.5 

1.4  X  1010 

5.  6  X  109 

2.5  X  1010 

8.  OX  109 

67 

7.5 

10 

8.4  X  109 

2.  7X  109 

1.7  X  1010 

4. OX  109 

130 

9.7 

- 

6.0  X  109 

- 

1 . 33  X  1010 

- 

200 

10 

20 

5.  7  X  109 

3.  3x  109 

1.3  X  109 

6.2  X  109 

210 

13.9 

- 

3.8  x  109 

- 

O 

o 

- 

370 

20 

30 

2.4  x  )09 

9.  OX  108 

6.8  10 

2.6  X  109 

700 

30 

40 

9 

1.5  X  10 

4.8^  108 

4.2  >  109 

1.6  ->  109 

1 ,450 

40 

50 

1.02  X  109 

2.  6  X  108 

2.6  X  109 

7.  8>  108 

2, 400 

50 

75 

7.6  X  108 

3.  6  X  108 

1.82  ■  109 

9.4V  108 

3,500 

75 

100 

4.0  X  108 

2.0  >  108 

00 

o 

00 

00 

4.  2  X  108 

7,000 

100 

200 

2.0  X  108 

1 . 7  x  108 

4.6  X  108 

3.4  V  108 

11,500 

200 

300 

3.  1  X  I07 

2.  4  >  107 

1 . 26  >  108 

7.0  >  107 

40,000 

300 

- 

7.0  X  106 

- 

5.6  X  107 

- 

77,000 
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TABLE  X 

PREDICTED  SOLAR  FLARE  ALPHA  PARTICLE  FLUX  OVER  THE  PERIOD  1967- 1973 


Energy  Range 
(MeV) 

2 

Alpha  Particle  Flux  (a/cm  •  year) 

Minimum  Shielding 
ta  Stap  Alpha  Particles 
af  Energy  E] 
(mg/cm2) 

Average 

Maximum 

Integral 

(E  >E,) 

Differential 

(E,  <E  <E2) 

Integral 

(E  >  E,) 

Differential 

(e,<e<e2) 

Ei 

E2 

1 

2 

5. 0  X  1010 

2.  2  X  IO10 

5.0  X  lo'° 

2. OX  IO10 

0.51 

2 

3 

2.  8  X  1010 

1 . 2  X  IO10 

3. 0  X  IO10 

9.  OX  109 

1.45 

3 

4 

1 . 6  X  1010 

5. 0  X  109 

2.  1  X  IO10 

5.  OX  109 

2.8 

4 

5 

1.1  x  io10 

4. 0  X  109 

1.6  X  IO10 

4. OX  109 

4.5 

5 

7.5 

7.  OX  109 

2.  8  X  109 

1 . 2  X  IO10 

3.  5  X  109 

6.75 

7.5 

10 

4.  2  X  109 

1.4X  109 

8.  5  X  109 

2.0  X  109 

13 

9.7 

- 

3. 0  X  109 

- 

6.  6  X  109 

- 

20 

10 

20 

2.  8  X  109 

1 . 6  X  109 

6.  5  X  109 

3.  1  X  109 

21 

13.9 

- 

1 . 9  X  109 

- 

5.0  X  109 

- 

36 

20 

30 

1.2X  109 

4. 5  X  108 

3.  4  X  109 

1.3X  109 

67 

30 

40 

7.  5  X  108 

2. 4  X  108 

2. 1  X  109 

8.0  X  108 

130 

40 

50 

5. 1  X  108 

1.3X  108 

1.3X  109 

3. 9  X  108 

210 

50 

75 

3.  8  X  108 

1 . 8  X  1 08 

9. 1  X  108 

4.  5  X  108 

310 

75 

100 

2. 0  X  108 

l.OX  io8 

4. 4  X  108 

2.  1  X  108 

625 

100 

200 

1.0X  io8 

7.  9  X  107 

2.  3  X  108 

2. 0  X  108 

1,050 

200 

300 

2. 1  X  107 

1.5X  107 

2.  9  X  107 

2.  3  X  107 

3,500 

300 

- 

6. OX  106 

- 

6. 3  X  106 

- 

7,000 
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V.  GALACTIC  COSMIC  RAYS 


A.  Composition  and  Spectrum 

The  partiele  flux  of  the  galaetic  cosmic  rays  is  small  eompared  to  that  of  the  trapped  radia¬ 
tion  and  of  solar  flare  partieles,  but  will  be  deseribed  here  for  the  sake  of  completeness.  Most 

of  the  data  on  the  charge  composition  and  speetra  of  the  various  components  of  the  galactic  eosmic 

113 

rays  have  been  taken  from  Webber's  recent  comprehensive  review.  Figure  27  shows  the  pri¬ 
mary  spectrum  of  the  total  galactic  radiation  over  the  rigidity  range  from  7  x  10*  to  10*  ^  V.*  *"*'* *^ 
At  the  low  rigidity  end  the  intensity  increases  by  about  a  factor  2  from  sunspot  maximum  to  sun¬ 
spot  minimum.  Over  a  wide  range  the  spectrum  may  be  expressed  by  a  power  law  as  a  function 
of  energy  or  rigidity. 

J(>E)  =  KEr  . 

9  15 

Over  the  entire  range  5  x  10  to  5  X  10  eV 

-1.60  ±0.03 

The  differential  primary  spectrum  obtained  by  differentiating  the  spectrum  of  Fig.  20  goes  through 
a  maximum  at  1  GY.  At  lower  energies  accurate  average  data  for  the  primary  radiation  are  dif¬ 
ficult  to  obtain  due  to  uncertainties  of  the  albedo  corrections  and  geomagnetic  cutoff  and  also 
due  to  the  large  effect  of  the  solar  modulation.  Such  data  have  been  obtained  recently  from  sat¬ 
ellite  measurements,  but  results  over  an  entire  sunspot  cycle  arc  not  yet  available. 

The  primary  galactic  cosmic  rays  are  largely  composed  of  protons  with  helium  nuclei  next 
in  abundance.  On  a  rigidity  basis  the  proton-helium  ratio  is  ~7.0.  The  integral  proton  spectrum 
may  be  obtained  directly  from  Fig.  27  by  decreasing  the  ordinate  by  the  18  percent  contribution 
of  the  heavier  nuclei.  The  more  recent  balloon  measurements  with  Cerenkov-seintillation  count¬ 
ers**^  shown  in  Fig.  28  exhibit  somewhat  higher  proton  intensities.  In  the  rigidity  range  from 
3  to  20  GY  the  proton  spectrum  may  be  expressed  as 

J(>P)  =  7000  P  *'^p/m^  •  sr  ■  sec 

where  P  is  in  gigavolts.  The  differential  proton  spectrum  at  low  rigidities  lias  been  obtained  by 

113 

Webber  from  observations  in  1963.  It  is  shown  in  Fig.  29  together  with  the  results  of  balloon 

and  satellite  measurements  carried  out  at  sunspot  minimum  by  Balasubrahmanyan*  *^  and  electron 

118 

measurements  by  L'Heureux  and  Meyer.  The  proton  spectrum  has  a  few  peaks  at  about  1  GV 
and  falls  off  sharply  at  lower  rigidities. 

The  integral  helium  spectrum  shown  in  Figs.  27  and  28  is  similar  to  the  proton  spectrum  and 
may  be  given  the  form 

.I(>P)  ~  900  P  *'S  helium  nuelei/m2  •  sr  •  sec 

over  the  rigidity  range  from  5  to  20  GV.  The  differential  helium  spectrum  peaks  al  1.3^  GV. 

Between  2.5  and  12  GY  the  differential  rigidity  spectra  of  helium  and  all  heavier  nuclei  may  be 

represented  by  power  laws  with  exponent  —2.2.  At  higher  rigidities  the  exponent  increases  to 

119 

—  2.5  or  -2.6  above  100  GV.  The  abundance  of  the  heavier  nuclei  is  listed  in  Table  \  .  Sec.  I\ 
Their  integral  spectra  are  similar  to  those  of  the  protons  and  of  helium,  their  differential  spectra 
peak  at  somewhat  higher  rigidities.  The  d/p  ratio  is  less  than  6  percent  in  the  energy  i  ange 
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25  —  80  MeV/nucleon;  about  30  to  40  percent  of  the  helium  nuclei  consist  of  He  on  a  rigidity 

114 

basis.  The  present  state  of  knowledge  on  primary  cosmic  ray  electrons  and  positrons  has 

115  114 

been  summarized  by  Meyer  and  VVaddmgton.  Between  500  Me\  and  3.5  GeY.  L'Houreux 
118 

and  Meyer  obtained  a  differential  spectrum  of  the  form 

dJ  ....  - 1.6/2 

=  11  x  E  e/m  .  sec  .  sr  .  GeV 

and  by  inference  an  integral  spectrum  with  exponent  —0.6.  An  extrapolation  of  this  spectrum  to 
lower  energies  is  in  agreement  with  previous  experimental  results.  Webber  and  Chotkowski* 
obtained  a  differential  spectrum 

-  30  x  E  * ' *  e/m2  *  sec  ■  sr  •  GeV 

Between  16  and  50  GeV,  Daniel  and  Stephens*2*  obtained  an  integral  spectrum  exponent  -  1.16. 
About  30  to  40  percent  of  the  particles  are  positrons.  The  fraction  of  positrons  seems  to 
decrease  with  energy. 

Gould  and  Burbidge*22  have  reviewed  the  experimental  measurements  on  high  energy  photons 

from  cosmic  sources.  Figure  30  shows  the  spectrum  of  the  isotropic  cosmic  background  photon 

fluxes  extrapolated  from  a  few  observations  at  different  energies.  Figure  31  shows  the  differ- 

123 

ential  spectrum  taken  from  the  paper  by  Rocchia,  et  al.  In  addition  discrete  galactic  x-ray 
sources  produce  fluxes  on  the  order  of  10  to  10  ergs/cm  •  sec. 

The  integral  flux  of  galactic  protons,  helium  nuclei  and  electrons  at  sunspot  minimum  is 
shown  in  Table  XI  for  different  energies.  The  radiation  is  very  low  in  intensity  but  very 
penetrating. 

B.  Temporal  Variations 

The  galactic  cosmic  radiation  is  modulated  by  the  sun  resulting  in  three  major  (  lasses  of 
126,127 

time  variations 

(1)  Solar  Cycle  Variations:—  The  intensity  varies  in  anticorrelation  with 
the  sunspot  activity  lagging  by  about  9  months  during  increasing  solar 
activity  and  by  18  months  during  periods  of  decreasing  solar  activity. 

The  modulation  is  a  function  of  rigidity  and  takes  the  form  oxp(  K  P)  in 
the  rigidity  range  0.5  P  <  50  GY. 128  At  1  GY  the  cosmic  ray  intensity 
increases  by  a  factor  2  from  sunspot  maximum  to  sunspot  minimum. 

(2)  Forbush  Decreases:—  Forbush  decreases  in  the  cosmic  ray  intensity 
on  the  order  of  20  percent  >1  GV  usually  occur  in  association  with 
magnetic  storms  and  last  from  2  days  to  2  weeks.  Small  decreases 
are  correlated  to  the  27-day  cycle  of  magnetic  activity. 

(3)  Local  Time  Variations:—  Very  small  diurnal  variations  (<0.5  percent 

1  GV)  take  place  in  the  cosmic  ray  intensity  with  an  average  maximum 
between  1300  and  1500  hours. 

C.  Geomagnetic  Cutoff 

The  geomagnetic  cutoff  rigidity  defined  in  Eq.  (15),  Sec.  IV.  is  a  more  meaningful  param¬ 
eter  for  galactic  cosmic  rays  during  sunspot  minimum  than  for  solar  particles.  Equation  (15) 
would  indicate  a  cutoff  rigidity  of  340  MY  or  a  cutoff  energy  of  60  MeY  at  L  t>.6.  This  is  below 
the  peak  energy  of  the  primary  r  adiation,  so  that  almost  the  entire  primary  cosmic'  ray  intensity 
will  reach  the  equatorial  synchronous  orbit.  Even  under  these  conditions  the  cutoff  rigidity  is 
subject  to  the  following  uncertainties: 
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TABLE  XII 

CHARGED  PARTICLE  RADIATION 

Particles 

Energy  Range 

Total  Dase  per  Year  in 
Equatorial  Synchronous  Orbit 

(particles/cm^) 

Trapped  Electrons 

<5  eV 

io17 

>200  eV 

1017 

>10  keV 

2.  1  X  1015 

>100  keV 

i.i  x  io'5 

>1  MeV 

1 . 5  X  1013 

Trapped  Protans 

<5  eV 

io17 

>200  eV 

1°17  ( 

>1.8  keV 

6  Xio'5 

>100  keV 

14 

3.  8  X  10 

>1  MeV 

5. 4  X  1010 

Salar  Pratans  (maximum) 

>1  MeV 

i  x  io11 

>10  MeV 

1.3X  lo'° 

>100  MeV 

4.  6  X  108 

Solar  Alpha  Particles 

>1  MeV 

5  X  1010 

(maximum) 

>10  MeV 

6.  5  X  109 

>100  MeV 

2.  3  X  108 

Galactic  Cosmic  Rays 

>100  MeV 

1. 7X  108 

(maximum) 

TABLE  XIII 

EFFECT  OF  SHIELDING  ON  CHARGED  PARTICLE  RADIATION 
IN  EQUATORIAL  SYNCHRONOUS  ORBIT 


Total  Dase  per  Year* 

(particl 

es/cm  ) 

6-mil  Si02 

7-mil  Kavar 

„  „  .14 

„  14 

Electrons 

3.9X  10 

1.3X  10 

Pratans 

2  X  1010 

8  X  109 

Alpha  Particles 

2.  5  X  109 

1  X  io9 

*  Omnidirectional  flux  with  infinite  shielding  aver  2tt  steradians. 
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(1)  No  experimental  determinations  of  the  cutoff  rigidity  in  synchronous 
orbit  have  been  published.  L,  =  6.6  may  not  be  a  good  parameter. 

129 

(2)  The  geomagnetic  cutoff  is  not  sharp. 

(3)  At  low  altitudes  Eq.  (15)  should  be  replaced  by 

P  =  KLy 
c 

where  K  varies  from  IS. 02  to  16.69  and  y  is  close  to —2. 

VI.  SUMMARY 

Table  XII  summarizes  the  charged  particle  radiation  in  the  equatorial  synchronous  orbit. 

It  may  be  noted  that  there  are  large  doses  of  very  low  energy  radiation  requiring  only  minimal 
shielding,  whereas  the  more  penetrating  radiation  is  far  less  intense.  While  it  is  beyond  the 
scope  of  this  report  to  consider  the  effect  of  different  shielding  configurations  on  the  incident 
charged  particle  spectra  (see  Refs.  131,  132),  two  situations  are  of  special  interest  to  communi 
cation  satellites  in  synchronous  orbit: 

(a)  6-mil  SiO^,  the  minimum  shielding  applied  to  the  silicon  solar  cells. 

(b)  7-mil  Kovar  or  nickel,  the  minimum  shielding  provided  by  the  can  or 
cover  of  a  hermetically  sealed  simiconductor  device. 

Table  XIII  indicates  that  even  these  two  minimum  shielding  conditions  remove  all  charged 
particles  except  trapped  electrons  as  well  as  some  solar  flare  protons  and  alpha  particles.  Ad 
ditional  shielding  information  is  provided  in  the  tables  in  the  text.  It  should  be  noted  that  the 
shielding  itself  creates  bremsstrahlung  which  may  contribute  a  significant  amount  of  ionizing 
radiation  when  most  of  the  primary  radiation  has  been  cut  off. 
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Fig.  1.  Integral  spectrum  af  trapped  electrons  in  synchronous  orbit  overoged  over  locol 
time  (taken  from  J.  |.  Vette  ond  A.  B.  Lucero,  Ref.  15). 


Fig.  2.  Comporison  of  energy  spectro 
of  tropped  electrons  in  synchronous  orbit 
(taken  from  W.  D.  Brown,  Ref.  19). 
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Fig.  3.  LES-4  electron  spectra. 
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Fig.  4.  Omnidirectional  electron  flux  integral 
spectrum  (based  on  calculations  by  J.l.  Vette 
on  OGO-A  spectrometer  doto). 
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Fig.  5.  Variotion  af  electron  flux  with  L. 


Fig.  6.  Variation  af  L  with  longitude  at 
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Fig.  7.  Temporal  variations  seen  by  LES-4  from  21  De¬ 
cember  1965  ta  8  July  1966. 


Fig.  8.  Composite  local  time  variation  (taken 
from  Ref.  15). 
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Fig.  9.  Pitch  angle  distribution  af  electrons  in 
synchronous  arbit  (taken  from  T.W.  Lezniak, 
et  ol . ,  Ref.  6). 


Fig.  11.  Omnidirectionol  proton  flux  integrol  "e 
spectrum  (based  an  calculations  by  J.  I.  Vette 
on  data  by  L.  R.  Dovis  and  J.  M.  Will  iamson,  * 
Ref.  52).  i 


Fig.  10.  Integral  spectrum  of  tropped  protons  in 
synchronous  orbit  (data  abtoined  from  J.  H.  King, 
Ref.  55). 
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Fig.  12.  Effect  of  solar  eventaf  5  October  1962  an 
trapped  proton  distribution  (token  from  L.  R.  Davis 
and  J.M.  Williamson,  Ref.  52,  Fig.  11). 
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Fig.  13.  Temporal  variation  of  trapped  proton 
intensities  meosured  on  Explorers  14  and  15 
(taken  from  L.  R.  Davis,  Ref.  51,  Fig.  12). 


Fig.  14.  Equatariol  pitch-ongle  distribution  of 
trapped  pratans  at  L  =  6  (taken  from  L.  R.  Davis, 
Ref.  51,  Fig.  8). 
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Fig.  15.  Differential  energy  spectrum  of  low  energy  trapped  protons  (based  on  data  from 
L.  A.  Frank,  Ref.  39  and  L.  R.  Davis,  Ref.  51). 
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Fig.  16.  Integrated  P0  distribution  for  solar- 
flare  events  (taken  from  W.  R.  Webber,  Ref.  67, 
Fig.  19). 


Fig.  17.  Proton-to-alpha-particle  ratio  as 
a  function  of  characteristic  rigidity  P0 
(taken  from  W.  R.  Webber,  Ref.  65,  Fig.  3). 
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Fig.  18.  Solor  x-roy  spectrum  (token  from  H.  H.  Malitsan,  Ref.  88). 
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Fig.  19.  Integrated  annual  intensity 
af  salar  particles  as  a  function  of 
average  yearly  sunspot  number  (based 
on  data  fram  W.R.  Webber,  Ref.  65). 
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PREDICTED  AND  OBSERVED  SUNSPOT  NUMBERS 


Fig.  20.  Predicted  ond  observed  sunspot  numbers  (token  from  Ref.  97). 
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Penetration  of  flare  particles  into  the  magnetosphere. 
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Fig.  22.  Integrol  onnual  solar  proton  intensity  vs  rigidity. 
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Fig.  23.  Differential  and  integral  pratan 
spectra  far  Bailey's  typical  solar  flare 
(taken  from  J.W.  Keller,  Ref.  110). 


Fig.  24.  Integrol  protan  flux  far  a  typical 
large  flare  (taken  fram  Ref.  111). 
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Fig.  26.  Predicted  integrated  salar  flare 
alpha  particle  spectra. 
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Fig.  27.  Integral  spectrum  af  galactic  cosmic 
rays  (total  primary  radiation  and  helium  nuclei 
from  W.R.  Webber,  Ref.  113,  electrons  fram 
C.J.  Waddingtan,  Ref.  114  and  P.  Meyer, 
Ref.  115). 


Fig.  28.  Integral  spectra  af  galactic  casmic  rays  (taken  fram  J.F.  Ormes 
and  W.  R.  Webber,  Ref.  1 16). 
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Fig.  29.  Differentiol  spectro  of  goloctic  cosmic  rays 
ot  low  rigidities  (protons  ond  helium  nuclei  1963  from 
W.  R.  Webber  ond  J.F.  Ormes,  Refs.  113,  116;  protons 
and  helium  nuclei  1965  from  V.  K.  Bolosubrohmonyon, 
et  aj.,  Ref.  117;  electrons  1964  from  J.  L'Heureux 
ond  P.  Meyer,  Ref.  118). 
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Fig.  30.  High  energy  cosmic  ray  photon 
spectrum  (taken  from  R.  J.  Gould  and 
G.R.  Burbidge,  Ref.  122). 
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Fig.  31.  Differential  cosmic  ray  photon  spectrum  (taken  from  R.  Rocchia,  et  al.,  Ref.  132). 
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